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disturbances  not  related  to  local  forcing  are  observed  to 
propagate  poleward.   At  6.1°N,  main  thermocline  displacements 
are  observed  to  propagate  westward  as  the  result  of  vertical 
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advection  of  mean  temperature  (colder  to  the  south)  by  sur- 
face Ekman  currents.   Eastern  boundary  disturbances  also 
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disturbances  observed  in  this  investigation  are  considered 
to  be  permanent  features  of  the  model  circulation.   Grid  size 
limitations  preclude  detailed  investigation  of  the  observed 
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I.   INTRODUCTION 

Recent  modeling  investigations  of  the  strength  and 
thermal  structure  of  Tropical  North  Pacific  Ocean  currents 
have  provided  evidence  as  to  the  existence  of  a  baroclinic 
response  to  annual  variations  in  wind  stress.   Meyers  (1978) 
has  observed  variations  in  the  depth  of  the  14°C  isotherm 
which,  at  6  N,  propagate  westward  at  nearly  the  speed  of 
free,  non-dispersive  Rossby  waves.   His  model  mechanisms  at 
6  N  are  Ekman  layer  divergence  (Ekman  pumping)  as  a  response 
to  local  wind  stress,  and  planetary,  geostrophic  divergence 
inherent  on  a  rotating  sphere.   At  10°N  Meyers  observed 
thermocline  depth  variations  which  were  nearly  in  phase 
across  the  Pacific  with  the  Ekman  pumping  velocity  calcu- 
lated from  the  local  wind  stress.   White  (1977)  found  an 
agreement  in  phase  only  at  150  W,  and  suggested  a  baroclinic 
response  to  annual  forcing  in  the  form  of  a  long  wave  which 
results  from  the  superposition  of  baroclinic  Rossby  waves 
emanating  from  the  eastern  boundary,  and  the  local  forced 
response.   This  resulting  long  wave  propagated  at  twice  the 
speed  of  free  Rossby  waves. 

The  purpose  of  this  study  was  to  investigate  further 
any  low  frequency  responses  to  wind  forcing  in  the  North 
Equatorial  Pacific  and  to  consider  whether  any  such  responses 
are  consistent  with  theory.   The  numerical  model  used  is  that 
of  Haney  (1974)  with  subsequent  improvements  involving 
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parameterization  of  convective  mixing  and  surface  wind  stir- 
ring (Haney  and  Davies,  1976)  and  nonlinear  eddy  viscosity 
(Haney  and  Wright,  1975).   This  model  is  based  on  the  hydro- 
static and  Boussinesq  approximations  and  is  a  10-level  prim- 
itive equation  model  in  a  closed  rectangular  basin  with  a 
uniform  depth  of  4  kilometers.   The  model  basin  extends  from 
the  equator  to  65  N  and  westward  for  90   of  longitude  from 
125  W.   Horizontal  resolution  is  approximately  3   of  longi- 
tude and  2°  of  latitude  with  33  uniformly  spaced  grid  points 
in  both  the  east-west  and  north-south  directions.   The  ten 
vertical  levels  are  placed  at  10,  32.5,  62.5,  102.5,  162.5, 
462.5,  900,  1700,  and  3100  meters  below  the  undisturbed  sea 
surface.   The  model  has  been  developed  and  improved  in  order 
to  study  large-scale  thermal  anomalies  observed  in  the 
Central  North  Pacific  Ocean  by  the  North  Pacific  Experiment 
(NORPAX). 

The  model  was  "spun  up"  for  a  number  of  years  of 
simulated  time  sufficient  to  allow  circulations  and  thermal 
patterns  to  establish  themselves  and  then  detailed  analysis 
of  the  transient  response  in  the  tropical  regions  was  made. 
The  wind  stress  prescribed  as  a  forcing  function  in  the  model 
from  the  equator  to  30  N  was  that  used  by  Meyers  (1978)  and 
described  in  detail  by  Wyrtki  and  Meyers  (1975).   North  of 
30°N  wind  stress  was  calculated  from  geostrophic  wind  data 
obtained  from  an  NCAR  data  tape  of  the  Northern  Hemisphere 
climatological  atlases  of  Jenne  et_  al_.  ,  (1969a,  b).   The 
major  region  of  interest  was  south  of  20  N.   This  location 
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was  where  Meyers  and  White  observed  wavelike  disturbances 
and  it  was  reasonable  to  expect  similar  waves  in  this  study 
By  analyzing  any  observed  waves  as  to  their  phase  relation- 
ship with  the  local  wind  stress  curl,  their  phase  speed, 
and  their  wavelength,  possible  mechanisms  for  their  genera- 
tion and  propagation  may  be  suggested  and  compared  with 
theory . 

By  better  understanding  how  the  oceanic  circulation 
changes  as  a  result  of  variant  wind  forcing,  more  compre- 
hensive studies  of  large-scale  air-sea  interaction  and 
climatic  changes  may  be  undertaken. 
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II.   DATA  SOURCES  AND  PROCESSING 

Data  from  two  different  sources,  one  covering  the  region 
north  of  30°N,  and  one  covering  the  region  south  of  30°N, 
were  combined  to  produce  a  realistic  wind  stress  field  which 
was  used  as  a  forcing  function  in  this  modeling  investiga- 
tion.  The  first  data  base  was  that  processed  and  explained 
in  detail  by  Wyrtki  and  Meyers  (1975).   This  monthly  aver- 
aged wind  stress  data,  sorted  into  2°  latitude  by  10°  longi- 
tude fields  encompassing  29  S  to  29  N  and  125  E  to  75  W, 
was  generously  provided  by  Dr.  Steve  Piascek  of  NORDA .   Only 
the  data  which  was  appropriate  for  use  by  the  model,  i.e. 
that  data  from  the  equator  northward  and  from  145  E  to  125  W, 
were  used  in  this  study.   Wyrtki  and  Meyers  calculated  these 
stress  data  for  each  observation  according  to  the  formula 

{xA,  x*}  =  pa  cd  |Va|  (ua,  va) 

where  x  is  the  wind  stress,  rectangular  component  positive 
to  the  east  (A)  and  north  (<$>),    p   is  the  constant  surface  air 

density  of  1.2  kg/m  ,  c,  is  the  drag  coefficient  assumed  con- 

_3 

stant  in  the  trade  wind  zone  at  1.5  x  10   ,  u   and  v   are 

'   a       a 

i  "*"  i       2      2  — 
the  surface  wind  components,  and  V   =  (u    +  v   )2,  the 

a       a      a 

magnitude  of  the  surface  wind.   A  quadratic  interpolation 
scheme  was  used  to  transfer  the  Wyrtki-Meyers  monthly  wind 
stress  component  arrays,  extending  from  the  equator  to  29°N, 
to  the  model  gridpoints. 
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The  second  source  of  data  was  an  NCAR  tape  of  monthly 
mean  climatological  surface  geostrophic  wind  component  data 
(Jenne  et_  al_.  ,  1969a,  b).   These  wind  components  had  been 
sorted  into  5°  latitude  by  5°  longitude  fields  encompassing 
the  model  latitudes  and  ranging  in  longitude  from  120  E  to 
80°W.   These  data  were  used  to  calculate  surface  wind  stress 

using  the  same  formula  as  above.   However,  an  atmospheric 

_3 
drag  coefficient  of  1.3  x  10    was  used  (Businger  et_  al.  , 

1971)  and  a  10%  reduction  in  magnitude  from  geostrophic  was 
assumed.   The  resulting  stress  components  for  the  model 
domain  north  of  30  N  were  interpolated  to  model  gridpoints. 

By  combining  the  stress  component  fields  computed  and/or 
interpolated  from  the  two  data  sources,  33  x  33  fields  en- 
compassing the  model's  geographical  domain  for  individual 

months  were  obtained.   The  Wyrtki-Meyers  stress  data,  in- 
cluded from  the  equator  to  29  N,  is  based  on  actual  wind 
observations  in  the  tropics,  while  the  NCAR  data,  from  30  N 
to  the  northern  boundary,  is  based  on  geostrophic  wind  cal- 
culations from  monthly  mean  sea-level  pressure  data.   It  is 
felt  that  the  resulting  stress  field  represents  the  best 
total  stress  pattern  that  is  presently  available  from  monthly 
mean  data  over  the  North  Pacific  Ocean. 

Although  the  major  feature  of  the  two  stress  fields 
matched  quite  well  at  the  common  latitude  of  30°N,  a  smooth- 
ing technique  was  employed  to  remove  small  scale  irregular- 
ities and  improve  horizontal  consistency.   A  weighted  three- 
point  meridional  average  given  by 
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T*i,rTi,j  +  -25(Ti,j+i  -2Ti,j  +  Ti,j-i>  wj 

was  selected  with  the  weight,  W, ,  a  Gaussian  function  of  the 
latitude  (i  is  the  longitudinal  gridpoint  index  and  j  is  the 

latitudinal  index).   W.  is  detailed  and  illustrated  in 

J 

Figure  1.   At  the  interface  between  the  two  data  sources, 
30  N,  W.  approaches  one,  and  at  the  latitudinal  boundaries 
of  the  model  it  approaches  zero.   Essentially,  then,  at  the 
interface  the  smoothing  process  was  a  \,    \,    \   meridional 
filter,  and  near  the  boundaries  was  negligible.   The  merid- 
ional smoothing  was  performed  twice  with  boundary  values 
being  unchanged.   Additionally,  a  i ,  J,  f  monthly  time 
filter  was  used  once  over  the  entire  fields  and  is  given  by 

t.  .  ,  =.25(t.  .  ,  ,  +  2t .  .  ,  +  t.    ,  ...  ) , 
i,J,k       i,j,k-l      i,j,k     i,j,k+l'' 

where  k  indicates  the  month. 

The  ocean  circulation  model  (Haney  et_  al_.  ,  1978)  was 
integrated  over  a  total  of  nine  complete  annual  cycles  start- 
ing from  the  results  of  the  previous  model.   The  only  change 
from  the  previous  model  is  the  introduction  of  the  new  wind 
stress  fields  as  described  above.   After  nine  years  of  inte- 
gration it  was  clear  that  the  upper  layers  of  the  ocean 
model,  i.e.  those  most  strongly  affected  by  the  annual 
variation  in  wind  forcing,  were  essentially  repeating  them- 
selves one  year  after  another.   The  results  analyzed  below 
are  those  from  the  eighth  and  ninth  years  of  the  nine  year 
integration.   Primary  attention  is  devoted  to  the  annual 
variation  of  model  variables,  less  to  the  annual  mean. 
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III.   OBSERVATIONS  AND  RESULTS 

The  trade  wind  field  over  the  Pacific  Ocean  is  one  of 
the  largest  and  most  consistent  wind  fields  on  earth  although 
little  is  actually  known  about  its  time  and  spacial  fluctua- 
tions (Wyrtki  and  Meyers,  1976).   Wind-driven  ocean  circula- 
tion studies  are  largely  dependent  on  wind  stress  field 
fluctuations  since  the  curl  of  the  wind  stress,  and  thus 
Ekman  surface  layer  divergence  or  pumping,  is  a  very  crucial 
forcing  parameter.   Accordingly,  maps  of  the  surface  wind 
stress  and  of  its  curl  were  analyzed  for  significant  fluc- 
tuations which  may  be  related  to  observable  oceanic 
disturbances . 

Figure  2  is  a  depiction  of  the  annual  mean  wind  stress 

vectors  over  the  model  domain.   The  region  of  primary  inter- 
est in  this  investigation  is  that  south  of  20  N.   As  shown, 
the  Intertropical  Convergence  Zone  (ITCZ)  has  an  annual 
mean  position  at  about  8  N.   Wyrtki  and  Meyers  (1975)  have 
shown  that  the  curl  of  the  wind  stress  and  thus  the  Ekman 
pumping  velocity  is  relatively  insensitive  to  changes  in  wind 
strength  and  depend  mostly  on  the  locations  of  strong  merid- 
ional shear,  which  show  latitudinal  variation  throughout  the 
year.   In  the  northern  hemisphere  winter  and  spring,  the 
northeast  trades  are  at  their  strongest  and  the  southeast 

trades  at  their  weakest.   The  ITCZ  lies  about  3   south  of 
its  annual  mean  position  and  as  seen  in  Figure  3,  a  map  of 
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wind  stress  curl  for  a  typical  spring  month,  regions  of 
strong  positive  wind  stress  curl  appear  between  5°N  and  10°N. 
Two  zones  of  particularly  strong  positive  curl  are  apparent, 
one  between  155°E  and  175°E  and  another  between  155  W  and 
140°W.   In  the  northern  hemisphere  summer  and  fall,  the 
ITCZ  shifts  northward  to  about  13°N  and  the  zone  of  strongest 
cyclonic  curl  is  located  between  10°N  and  15°N  as  seen  in 
Figure  4. 

The  strongest  meridional  shear  in  the  surface  wind  stress 
is  observed  in  a  band  slightly  north  of  the  ITCZ  which  passes 
over  the  annual  mean  position  of  the  ITCZ  twice  per  year. 
Therefore  the  amplitude  of  the  first  annual  variation  of  the 
shear  is  largest  at  the  latitudinal  extremes  of  the  ITCZ  and 
small  at  the  annual  mean  position.   Meyers  (1978)  computed 
the  amplitude  of  the  first  annual  harmonic  of  the  Ekman  pump- 
ing velocity  and  the  depth  of  the  14°C  isotherm  in  the  North 
Pacific  and  found  that  both  reached  maxima  at  6  N  and  at  10  N. 
Not  surprisingly,  these  latitudes  coincided  with  the  zones  of 
strongest  cyclonic  wind  stress  curl.   As  mentioned  previously, 
at  10  N  the  Ekman  pumping  velocity  and  the  depth  of  the  14°C 
isotherm  fluctuated  in  phase  across  the  Pacific  while  at  6°N 
the  temperature  fluctuations  propagated  westward  from  a 
region  of  resonant  forcing  in  the  Eastern  Pacific.   In  order 

to  analyze  results  which  could  best  be  compared  with  those 
of  Meyers  (1978),  model  latitudes  of  6.1°N  and  12.2°N  were 
chosen  for  study  in  this  investigation. 
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a.   Results  at  12.2°N 

Meyers'  (1978)  result  at  10  N  was  an  oceanic  tempera- 
ture response  which  was  nearly  in  phase  (a  two  month  zonal 
lag)  with  the  local  Ekman  pumping  velocity  across  the  entire 
Pacific  Ocean.   White  (1977)  theorized  that  the  annual  fluc- 
tuation at  any  interior  ocean  location  from  10  N  to  20  N 
depends  on  the  fluctuating  wind  field  over  the  entire  east- 
ward ocean  with  the  phase  speed  and  amplitude  strongly  related 
to  distance  from  the  eastern  boundary.   In  a  bounded  ocean, 
the  presence  of  the  eastern  boundary  can  significantly  alter 
the  phase  distribution  of  the  oceanic  response  to  Ekman  pump- 
ing by  acting  as  a  concentrated  source  of  baroclinic 
planetary  waves  which  propagate  boundary  responses  westward 
into  the  interior,  drastically  altering  the  phase  relation 
predicted  by  local  forcing  only  (White,  1977). 

Model  level  5,  corresponding  to  a  depth  of  162.5 
meters,  was  chosen  for  detailed  study  since  it  is  the  model 
level  which  most  closely  corresponds  with  the  observed  depth 
of  the  14°C  isotherm  (Meyers,  1978).   This  level  is  above  the 
intermediate  layers  where  the  vertical  temperature  gradient 
is  weak  and  below  the  level  of  the  maximum  density  gradient, 
whose  stability  shields  it  from  temperature  fluctuations 
associated  with  sea  surface  heat  exchanges.  Time-longitude 
maps  of  departures  from  the  annual  mean  (denoted  with  a  prime) 
of  the  zonal  (u')  and  meridional  (vf)  flow,  vertical  velocity 
(w'),  and  temperature  (_T '  )  were  developed  to  show  any  zonal 
propagation  of  wavelike  disturbances.   Figures  5  and  6  are 
such  maps  of  u'  and  v'  respectively  at  level  5  across  12.2°N 
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over  the  last  two  years  of  the  model  run.   These  figures  in- 
dicate westward  propagating  annual  signals  which  emanate  from 
near  or  at  the  eastern  boundary  with  the  zonal  component  hav- 
ing a  slightly  faster  phase  speed  and  a  higher  magnitude. 
Figure  7  is  a  similar  depiction  of  the  disturbance  vertical 
velocity,  w' .   The  major  features  of  the  open  ocean  vertical 
velocity  field  in  Figure  7  develop  near  the  eastern  boundary 
and  seem  to  propagate  westward  as  far  as  the  mid  ocean  region 
at  which  point  further  westward  propagation  with  increasing 
time  is  not  observed.  White  (1977)  observed  that  due  to  the 
inverse  latitudinal  dependence  in  zonal  phase  speed,  baro- 

clinic  long  waves  in  the  eastern  ocean  are  directed  westward 
while  in  the  western  ocean  they  appear  to  be  directed  more 
northward  due  to  refraction  effects.   Regardless,  the  overall 
pattern  observed  in  Figure  7  does  not  correlate  with  those 
of  the  zonal  and  meridional  flow  components  shown  in  Figures 
5  and  6  respectively.   Vertical  velocity  is  generally  related 
to  the  divergent  part  of  the  flow;  the  propagating  features 
of  Figures  5  and  6  are  most  likely  due  to  the  rotational  part 
of  flow  and  as  such  have  little  impact  on  the  vertical 
velocity.   It  is  noted  that  disturbances  in  Figures  5  through 
7  appear  to  intensify  along  the  western  boundary  of  the  model 
ocean.   This  is  to  be  expected  since  westward  propagating 
features  tend  to  congregate  at  or  near  western  boundaries 
in  models  such  as  that  used  in  this  study. 

Figure  7  also  shows  local  disturbances  of  great  inten- 
sity along  the  eastern  boundary.   However,  these  disturbances 
do  not  appear  to  propagate  westward  to  any  significant  extent, 
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They  also  do  not  appear  to  be  related  to  the  longshore  (merid- 
ional) component  of  wind  stress  nor  to  local  Ekman  pumping. 
In  order  to  investigate  any  meridional  propagation  of  these 
boundary  features,  a  time-latitude  map  of  w'  along  the  eastern 
boundary  was  prepared  and  is  shown  in  Figure  8.   Wyrtki  (1975) 
has  surmised  that  the  annual  baroclinic  response  at  an  eastern 
boundary  away  from  the  equator  is  in  part  composed  of  a  per- 
turbation which  propagates  into  the  region  from  the  equatorial 
region  in  the  form  of  an  internal  Kelvin  wave.   Those  coastal 
disturbances  at  12.2  N  in  Figure  8  have  a  computed  (northward) 
phase  speed  of  11.7  cm/s  and  the  major  features  have  an  annual 
period.   It  is  interesting  to  note  that  Enfield  and  Allen 
(1978)  have  recently  observed  spatially  coherent  disturbances 
in  monthly  mean  anomalies  of  temperature  and  sea  level  along 

the  Pacific  coast  of  North  and  South  America.   In  the  tropical 
latitudes,  these  observed  fluctuations  are  also  uncorrelated 
with  local  forcing  mechanisms  and  propagate  poleward  at  phase 
speeds  comparable  to  those  obtained  in  the  model  (Allen, 
personal  communication).   Because  of  the  coarse  grid  in  the 
present  model  however,  no  further  comparison  between  model 
results  and  coastal  observations  is  justified. 

Another  important  result  which  can  be  seen  from  Figures 
5  through  8,  and  others  below,  is  that  the  two  consecutive 
model  years  are  virtually  identical.   This  means  that  the 
model  is  in  a  statistical  equilibrium  with  the  forcing  and 
that  the  oscillations  being  examined  are  not  transitory  but 
are  permanent  features  of  the  model  circulation. 
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Figure  9  is  a  time-longitude  map  of  temperature  dis- 
turbances which  is  very  consistent  with  the  vertical  velocity 
pattern  given  in  Figure  7.   A  two  month  phase  lag  is  apparent 
with  generally  warm  water  sinking  during  the  first  half  of 
the  year  and  cold  water  rising  during  the  second  half  - 
clearly  and  indirect,  thermally  forced  disturbance.   Calcula- 
tions based  on  the  annual  positive  disturbance  in  the  east- 
central  ocean  in  Figure  9  revealed  a  zonal  wavelength  of 
about  10,000  kilometers  and  a  zonal  phase  speed  of  approxi- 
mately -42  cm/s.   These  results  are  in  qualitative  agreement 
with  White  (1977)  who  theorized  a  wavelength  of  10,000  kilo- 
meters and  a  phase  speed  of  -32  cm/s  and  observed  a  wavelength 
of  16,000  kilometers  and  a  phase  speed  of  -46  cm/s  at  12  N. 

The  Ekman  pumping  velocity  was  calculated  from  the 
following  formula  as  described  in  detail  by  DeWitt  and 
Leetmaa  (1978) 

0 

W  ,  =  /    (u   +  v  )  dz. 
-h   J  x     y 

-h 
A  map  of  the  disturbance  Ekman  pumping  velocity  is  presented 
in  Figure  10.   This  pattern  corresponds  with  the  gross  fea- 
tures of  the  vertical  velocity  field  of  Figure  7  although  no 
propagation  is  noted  in  the  eastern  ocean.   This  is  in  con- 
trast with  the  results  of  both  White  (1977)  and  DeWitt  and 
Leetmaa  (1978)  who  noticed  a  tendency  for  the  Ekman  pumping 
disturbances  to  propagate  westward.   Magnitude  differences 
between  corresponding  features  in  Figures  7  and  10  are  perhaps 
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due  to  the  fact  that  Figure  10  represents  the  depth  of  the 
bottom  of  the  Ekman  layer  which  need  not  be  coincident  with 
level  5  of  the  model.   Additionally,  different  processes 
are  involved  in  vertical  velocity  and  Ekman  pumping.   The 
theoretical  Ekman  depth  was  calculated  according  to  the 
formula 

2K 

De  =  ^-r1)2 

(DeWitt  and  Leetmaa,  1978)  where  K   is  the  vertical  eddy  vis- 

2 

cosity  coefficient  in  the  model  which  is  10  cm  /s,  and  f  is 

the  Coriolis  parameter.   At  12.2°N,  this  depth  was  calculated 
to  be  approximately  20  meters. 

In  order  to  describe  the  vertical  structure  of  the 
annual  wave  disturbances,  vertical  cross  sections  ( x  ,  z  )  of 
the  disturbance  quantities  were  prepared  across  12.2  N  on 
August  3  of  the  last  model  year.   This  date  was  selected  so 
as  to  coincide  with  the  large  positive  amplitude  disturbance 
noted  on  the  time-longitude  section  of  Ekman  pumping  in 
Figure  10.   Figure  11,  the  vertical  cross  section  of  u"  on 
August  3,  may  be  compared  with  the  time-longitude  map  of  u * 
in  Figure  5.   The  positive  wavelike  disturbance  near  the 
eastern  boundary  in  Figure  5  extends  from  the  surface  to  the 
bottom  with  greatest  magnitude  near  the  surface,  as  is  the 
case  with  the  positive  disturbance  in  the  western  ocean.  The 
negative  disturbance  in  the  mid-ocean  extends  from  th.e  surface 
to  about  1500  meters  with  greatest  magnitude  at  about  50  meters. 
The  result  that  the  largest  magnitudes  are  not  at  the  surface 
in  Figure  11  may  be  explained  by  the  dependence  of  the  zonal 
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flow  on  temperature  and  pressure  gradient  structure  and  not 
so  much  on  the  surface  winds.   Generally,  then,  the  zonal 
disturbances  extend  quite  deep  at  their  formative  stage  and 
do  not  significantly  deepen  or  shoal  as  they  propagate  west- 
ward. 

Figure  12  shows  the  corresponding  cross  section  of 
the  meridional  velocity  perturbation  which  may  be  compared 
with  the  time-longitude  section  of  Figure  6.   The  positive 
disturbance  in  the  mid-western  ocean  in  Figure  6  on  August  3 
extends  vertically  to  within  about  50  meters  of  the  surface. 
The  negative  disturbance  in  the  eastern  ocean  extends  from 
about  1500  meters  to  the  surface  where  it  expands  to  cover 
the  entire  model  domain  across  12.2  N.   The  negative  distur- 
bance is  of  greatest  magnitude  near  the  surface  at  all  longi- 
tudes while  the  positive  disturbance  is  very  small  whenever 
it  appears.   Apparently  the  meridional  flow  component  is 
more  closely  related  to  the  surface  winds  than  is  the  zonal 
component . 

Figure  13  depicts  the  vertical-zonal  cross  section 
of  vertical  velocity  disturbance  across  12.2  N  on  August  3. 
The  positive  disturbance  in  Figure  7  encompasses  most  of  the 
model  domain  except  for  the  surface  layer  and  the  extreme 
eastern  boundary.   Largest  magnitudes  are  centered  at  depths 
of  about  1000  meters.   The  eastern  boundary  negative  distur- 
bance extends  from  about  200  meters  to  the  surface. 

Corresponding  temperature  disturbances  are  shown  in 
Figure  14.   Warm  conditions  are  noted  in  all  areas  on  August 
3  except  the  eastern  region  at  levels  below  about  75  meters. 
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Wunsch  (1977)  has  noted  that  if  wind  forcing  is  as- 
sumed to  be  over  an  annual  period  then  there  is  no  reason  to 
isolate  the  problem  of  vertical  propagation  from  that  of 
horizontal.   He  contends  that  internal  gravity  waves  will 
determine  the  vertical  modal  structure  and  that  in  the 
presence  of  a  non-dissipat ive  and  flat  bottom  energy  may 
travel  from  the  surface  to  the  bottom  and  then  be  reflected 
back  to  the  surface.   Any  vertical  propagation  of  wavelike 
disturbances  should  be  apparent  in  time-depth  maps  of  the 
perturbation  quantities.   These  depictions  were  developed 
and  are  shown  in  Figures  15  through  18.   These  maps  were  made 
at  an  arbitrarily  selected  interior  ocean  longitude  of  162  W. 
No  vertical  propagation  is  seen  in  the  zonal  and  meridional 
velocity  components,  Figures  15  and  16  respectively.   Again 
it  can  be  seen  that  u'  penetrates  deeper  than  v'  and  is 
therefore  most  likely  due  to  the  north-south  meandering  of 
the  major  equatorial  current  systems  which  are  primarily 
zonal.   Figure  17  reveals  that  the  vertical  velocity  distur- 
bances also  do  not  appear  to  propagate  vertically  but  they 
attain  their  greatest  magnitudes  at  fairly  substantial 
depths.   These  deep  centers  apparently  are  the  cause  of  the 
corresponding  deep  temperature  disturbances  (Figure  18) 
which  lag  the  vertical  velocity  by  about  2 \   months.   The 
temperature  fluctuations  above  150  meters  are  due  to  the 
seasonal  cycle  of  surface  heat  flux  modified  slightly  by  a 
fluctuating  meridional  temperature  advection  in  the  upper  50 
meters.   Again  notice  how  the  two  model  years  are  nearly 
identical  -  thus  in  equilibrium  with  the  forcing. 
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The  model  results  at  12.2  N  are  summarized  as  follows 
The  dominant  thermal  response  at  the  depth  of  the  main  thermo- 
cline  occurs  in  phase  across  the  entire  longitudinal  extent 
of  the  ocean  and  is  attributed  to  large  scale  vertical  advec- 
tion  produced  by  wind-induced  Ekman  pumping.   This  is  in 
agreement  with  the  theoretical  analysis  and  observational 
results  of  Meyers  (1978).   In  addition,  there  is  a  pronounced 
fluctuation  in  the  zonal  current  in  the  open  ocean  which 
appears  to  be  associated  with  a  north-south  seasonal  meander- 
ing of  the  main  zonal  current  systems.   There  also  exists 
poleward  propagating  disturbances  confined  to  the  immediate 
vicinity  of  the  eastern  boundary  which  do  not  appear  to  be 
related  to  the  local  forcing.   Further  study  is  required  to 
determine  whether  these  disturbances  are  in  any  way  related 
to  the  sea  level  disturbances  observed  along  the  North 
American  coast  by  Enfield  and  Allen  (1978). 
b.   Results  at  6 . 1°N 

Variations  in  the  depth  of  the  14  C  isotherm  which 
propagated  westward  at  nearly  the  phase  speed  of  free  non- 
dispersive  Rossby  waves  were  observed  at  6  N  by  Meyers  (1978) 
He  theorized  that  the  waves  were  generated  by  the  zonal  and 
seasonal  variability  in  Ekman  layer  divergence  in  the 
eastern  Pacific  and  that  they  propagated  nearly  freely  in 
the  western  Pacific  where  he  found  the  forcing  was  weak.  The 
resulting  wave  lagged  the  Ekman  pumping  velocity  by  slightly 
more  than  one  month  at  the  position  of  maximum  forcing. 


Figures  19  and  20  are  time-longitude  maps  of  zonal 
and  meridional  flow  disturbances,  uf  and  vT  respectively, 
across  6.1  N  at  162.5  m  (model  level  5)  during  the  same  two 
model  years  analyzed  above  at  12.2  N.   These  figures  indicate 
rapidly  westward  propagating  annual  signals  which  appear  to 
emanate  from  the  eastern  boundary.   The  zonal  component  in 
Figure  19  appears  to  increase  gradually  in  magnitude  to  a 
maximum  in  the  region  slightly  west  of  the  mid-ocean  and  then 
decreases  from  there  to  the  western  boundary.   The  zonal  dis- 
turbances propagate  faster  than  the  meridional  and  faster 
at  6.1°N  than  at  12 . 2°N  (Figure  5).   Figure  21  is  a  time- 
longitude  map  of  vertical  velocity  disturbance  at  the  same 
latitude  and  depth.   As  was  noted  for  Figure  7  at  12.2  N, 
local  disturbances  of  great  magnitude  are  observed  along  the 
eastern  boundary  which  do  not  appear  to  propagate  westward. 
These  northerly  propagating  Kelvin-like  disturbances  (it  can 
be  seen  from  Figures  19  through  22  that  v' ,  T' ,  and  wr  are 
involved  in  the  disturbance  but  not  uf)  are  depicted  in 
Figure  8.   Near  the  eastern  boundary  in  Figure  21,  wavelike 
perturbations  form  and  propagate  completely  across  the  model 
domain,  in  contrast  with  those  at  12.2  N  (Figure  7)  which 
appear  to  be  nearly  in  phase  across  the  ocean.   The  distur- 
bances' amplitudes  slowly  vary  with  westward  propagation. 

The  temperature  disturbance  at  6.1  N  (Figure  22) 
has  essentially  the  same  pattern  as  the  vertical  velocity 
disturbance  but  with  a  two  month  phase  lag  which  indicates 
the  temperature  is  responding  to  the  vertical  velocity.   Cal- 
culations revealed  a  phase  speed  of  -32  cm/s,  a  zonal 
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wavelength  of  about  9900  kilometers,  and  an  annual  period. 
Meyers  (1978)  observed  a  corresponding  phase  speed  of  -56 
cm/s  with  a  wavelength  of  17,000  kilometers  at  6  N. 

A  time-longitude  map  of  Ekman  pumping  velocity  distur- 
bance across  6.1  N  is  shown  in  Figure  23.   The  theoretical 
Ekman  depth  as  calculated  above  for  12.2  N  was  found  to  be 
approximately  28  meters  at  6.1°N.   The  pattern  in  Figure  23 
correlates  well  with  the  actual  vertical  velocity  (Figure  21) 
which  in  turn  correlates  with  the  temperature  fields  except 
for  the  region  of  very  slow  westward  propagation  in  the  mid- 
ocean  area.   Propagation  is  very  rapid  in  both  the  eastern 
and  western  thirds  with  greatest  amplitudes  observed  at  the 
eastern  boundary  (125  W).   The  latter  observation  is  consis- 
tent with  Meyers'  (1978)  finding  that  the  maximum  displace- 
ment occurs  at  130  W.   The  region  of  strong  cyclonic  wind 
stress  curl  observed  at  this  latitude  in  Figure  3  was  at 
and  slightly  west  of  140°W. 

The  zonal-depth  structure  of  the  disturbance  patterns 
across  6.1  N  on  August  3  of  the  last  model  year  is  depicted 
in  Figures  24  through  27.   Figure  24  shows  that  the  positive 
u'  disturbance  of  Figure  19  is  very  broad  at  great  depths, 
narrow  at  mid  levels,  and  widens  again  in  the  surface  layers 
with  slightly  higher  magnitudes  observed  at  the  surface.  The 
corresponding  negative  disturbance  is  widest  and  has  its 
greatest  magnitude  at  about  75  meters.   As  was  discussed  pre- 
viously with  regard  to  Figure  11,  there  is  a  deep  zonal 
current  fluctuation  which  seems  to  be  related  to  the  north- 
south  displacement  of  the  main  current  systems  which  are 
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zonal.   The  meridional  velocity  disturbance  across  6.1  N  is 
presented  in  Figure  25.   A  relatively  weak  flow  is  predomi- 
nant across  the  model  domain  except  in  the  surface  layers 
where,  as  in  Figure  12,  there  is  a  strong  response  due  to 
the  dependence  on  surface  winds.   Positive  disturbances  shown 
in  Figure  25  are  of  much  smaller  magnitude. 

The  zonal-depth  map  of  w'  in  Figure  26  shows  that 
with  the  exception  of  a  narrow  semi-vertical  band  in  the 
western  ocean  and  the  Kelvin-like  disturbances  on  the  eastern 
boundary,  the  model  domain  is  predominantly  occupied  by  sink- 
ing (w'  <  0)  disturbances  at  this  time.   The  shallow  Kelvin 
disturbance  extends  deeper  at  6.1  N  than  it  did  at  12.2  N 
(Figure  13)  and  a  phase  change  is  apparent  between  the  two 
latitudes.   Figure  27  shows  that  the  model  domain  is  colder 
than  normal  except  for  the  eastern  ocean  surface  layer  and 
a  mid-level  warm  disturbance  in  the  western  ocean. 

Time-depth  maps  of  u'  and  v'  at  162  ¥  are  shown  in 
Figures  28  and  29.   Disturbances  of  greatest  magnitudes  occur 

at  the  surface  with  those  of  the  zonal  component  extending 
quite  deep.   On  comparing  these  figures  with  their  counter- 
parts at  12.2°N,  Figures  15  and  16,  the  differences  in  dis- 
turbance magnitudes  between  the  two  latitudes  is  obvious;  the 
deep  response  is  a  little  stronger  at  6 . 1°N  but  the  surface 
Ekman  flow  is  nearly  an  order  of  magnitude  higher  at  6.1  N. 
Plots  of  the  disturbances  of  the  zonal  and  meridional  com- 
ponents of  wind  stress  as  functions  of  time  at  162  W  are 
shown  in  Figure  30.   The  large  magnitude  surface  layer  distur- 
bances noted  in  Figures  28  and  29  are  in  phase  with  these 
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stress  disturbances  indicating  that  the  wind  stress  is  most 
likely  the  dominant  forcing  process  of  ocean  (Ekman)  flow 
at  the  surface  and  immediately  below.   The  time-depth  map 
of  vertical  velocity  disturbance  across  6.1  N  as  presented  in 
Figure  31  is  similar  to  that  of  12.2°N  (Figure  17)  in  that 
greatest  magnitudes  are  reached  at  deep  levels.   However,  in 
contrast  to  12.2  N,  there  is  a  phase  shift  at  intermediate 
levels  at  6.1  N.   The  phase  shift  probably  represents  a 
transition  region  between  the  upper  500  meters  or  so,  which 
is  responding  to  the  wind,  and  the  remainder  of  the  ocean 
which  is  reacting  out  of  phase.   Figure  32  shows  the  strong 
response  of  the  temperature  field  in  the  upper  500  meters  to 
the  vertical  motion  with  a  two  month  lag.   At  the  same  time, 
surface  disturbances  due  to  seasonal  wind  changes  occur  but 
they  are  confined  to  the  upper  50  meters.   At  6 . 1°N  vr  advects 
cold  water  northward  from  the  equatorial  surface  region  during 
the  first  half  of  the  year  and  warm  water  southward  during  the 
second  half  of  the  year. 

The  model  results  at  6.1  N  are  summarized  as  follows. 
The  dominant  thermal  response  in  the  vicinity  of  the  main 
thermocline  (level  5,  or  162.5  meters  in  the  model)  is  a  west- 
ward propagating  disturbance  (Figure  22)  produced  by  vertical 
advection  (w'  in  Figure  21)  induced  by  Ekman  pumping  (Figure 
23).   Figure  27  shows  that  the  temperature  response  seen  at 
162.5  meters  involves  the  entire  thermocline  and  is  actually 
strongest  at  about  100  meters.  In  addition  to  the  propagating 
thermal  response  in  the  main  thermocline  there  is  an  equally 
pronounced  annual  cycle  of  surface  temperature  CFigure  32) 
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which  is  produced  primarily  by  meridional  advection  of  mean 
temperature  (colder  to  the  south)  by  surface  Ekman  currents 
(Figure  29). 
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IV.   CONCLUSIONS 

Wavelike  features  simulated  in  the  tropical  regions  of 
an  ocean  circulation  model  which  are  seen  as  disturbances 
from  the  annual  mean  of  vertical  velocity  and  temperature 
have  been  described.   These  features  are  considered  to  be 
permanent  oscillations  of  the  model  circulation  and  result 
from  variability  in  wind  stress  forcing.   At  12.2  N,  a 
thermal  response  to  large  scale  vertical  advection  produced 
by  wind-induced  Ekman  pumping  or  surface  layer  divergence 
occurs  practically  in  phase  across  the  entire  longitudinal 
extent  of  the  model  ocean.   Underlying  features  of  the 
annual  signal  do  appear  to  propagate  westward.   These  results 
are  in  agreement  with  the  theory  and  observations  of  Meyers 
(1978).   Evidence  also  suggests  a  zonal  fluctuation  in  the 
interior  ocean  current  systems.   Along  the  eastern  boundary 
poleward  propagating  perturbations  of  the  vertical  velocity 

field  exist  which  may  be  related  to  sea  level  disturbances 
which  have  been  observed  by  Enfield  and  Allen  (1973). 

At  6.1  N,  thermocline  displacements  propagate  westward 
also  as  a  result  of  vertical  advection  induced  by  Ekman  pump- 
ing.  Disturbances  at  this  latitude  are  significantly  higher 
in  magnitude  than  those  observed  at  12.2  N.   Additionally, 
pronounced  sea  surface  temperature  disturbances  are  apparently 
produced  by  meridional  advection  of  mean  temperature  by 
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surface  Ekman  currents.   Eastern  boundary  perturbations  also 
exist  and  extend  deeper  at  6.1  N  than  at  the  higher  latitude 

Results  and  observations  in  this  study  were  generally 
consistent  with  accepted  theory.   It  would  be  beneficial  to 
consider  latitudes  other  than  6.1  N  and  12.2  N  in  future 
studies  in  order  to  more  completely  describe  the  tropical 
wavelike  disturbances  which  were  noted.   The  equatorial 
boundary  in  the  current  model  may,  however,  cause  unreason- 
able results  at  lower  latitudes  than  6.1  N.   The  model  is 
presently  not  suitable  for  use  for  further  study  of  the 
eastern  boundary  disturbances  because  of  grid  size  limita- 
tions.  However,  a  careful  study  of  these  features  using  a 
higher  resolution  model  having  a  realistic  coastline  con- 
figuration may  be  valuable  in  light  of  the  recent  findings 
of  Enfield  and  Allen  (.1978). 

The  present  model  has  little  or  no  synoptic  "noisert. 
It  is  recognized  that  this  feature  may  be  altered  if  finer 
resolution  or  synoptic  wind  forcing  were  introduced. 
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